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ABSTRACT 



Context. We report the first spectroscopic mapping of an atomic carbon line in an infrared dark cloud (IRDC). 

Aims. By observing the spatial distribution of the [Ci] emission in an IRDC, comparing it with the 13 CO emission and the known 
distribution of internal heating sources, we can quantify the role of internal and external UV irradiation in the production of atomic 
carbon. 

Methods. We used the 2 x 4 pixel SMART receiver of the KOSMA observatory on Gornergrat to map the [C i] 3 P v - 3 P line in the 

IRDC G48.65-0.29 and compared the resulting spectra with data from the BU-FCRAO 13 CO 1-0 Galactic Ring Survey. 

Results. The [C i]/ 13 CO effective beam temperature ratio falls at about 0.3 with local deviations by less than a factor two. All velocity 

components seen in 13 CO are also detected in [C ij. We find, however, significant differences in the morphology of the brightest regions 

seen in the two tracers. While 13 CO basically follows the column density distribution derived from the near-infrared (NIR) extinction 

and the submm continuum, the [C i] emission peaks at the locations of the two known NIR point sources. We find C i/CO abundance 

ratios between 0.07 and 0.13 matching the lower end of the range previously measured in star-forming regions. 

Conclusions. Evaluating the relative importance of the irradiation by embedded sources and by the Galactic interstellar UV field, we 

find that in G48.65-0.29 most [C i] emission can be attributed to externally illuminated surfaces. Embedded sources have a significant 

impact on the overall abundance distribution of atomic carbon as soon as they are found in an evolved state with noticeable NIR flux. 

Key words. ISM: clouds - ISM: structure 



1. Introduction 

Chemical models of interstellar gas predict that atomic carbon, 
Ci, occurs only in a transitional layer characterised by a moder- 
ate UV field which allows ionised carbon to r e comb i ne, but still 
leads to the dissociation of CO dLe Bourlotl |1993t lOka et al.L 
2004; Rollig et al., 2006). In the picture of a simple cloud ge- 
ometry with a well-defined interior and exterior of a cloud, [C i] 
is thus supposed to trace surfaces of molecular clouds exposed 
to UV radiation either from the external interstella r radiation 
field or from embedded stars (Ossenk opf et alj,[2007l) . Its abun- 
dance and emission is modelled in the framework of photon- 
dominated regions (PDRs). Molecular cloud observations show, 
however, that their structure is rather described by a fractal dis- 
tribution, where UV radiation can deepl y penetrate and surface s 
are found throughout the whole cloud dFalgarone et all l2004t) . 
leading to a strong correlation in the distribution of atomic car- 
bon and molecular tracers. 

Bright [C i] emission has been observed towards severa l 
star-forming regions (e.g.. iPlume et all [20001 ISun et all 12008). 
where the intense UV radiation from young stars creates many 
PDRs throughout the clumpy ambient molecular clouds. To un- 
derstand the origin of the [C i] emission and the Galactic abun- 
dance distribution of atomic carbon it is, however, essential to 
distinguish the mutual importance of PDRs within molecular 
clouds, created by star formation, and the PDRs at their outer 
surfaces p roduced by the int e rstella r radiation field illuminat- 
ing them dDraine & Bertoldil 119961) . iPineda & Benschl d2007l) 
detected [Ci] towards B68, only seeing the standard interstel- 



lar radiation field, but didn't find a PDR model that fitted all of 
the observed lines. 

To help resolving the question on the relative importance of 
inner and outer PDRs for the production of atomic carbon we 
decided to map the [C i] emission from an infrared dark cloud 
(IRDC), a s ite of potentially massive star formation i n a very 



early stage (Rathb orne et 



ally massivs 
iil l2006H s 



imon et all 120061) . IRDCs 



were detected as dark patches in front of the bright diffuse IR 
emission f rom the Mi l ky Way mapped by ISO and the MSX 
satellite dPerault et all Il996t lEgan et all 1 1998b . Their column 
densities exceed 10 23 crrT 2 , producing the strong dust absorp- 
tion that led to their discovery. IRDC masses range from 10 3 
to 10 4 M©, the gas is typica l ly col d (T < 25 K) and dense 
(n > 10 5 cm 3 , ISimon et all 120061) . Recent follow-up obser- 
vations have reve aled signs of active star (clust er) formation in 
numerous IRDCs dRathborne et alll2006ll2007t) . 

We have observed the IRDC G48.65-0.29 (ISimon et all 
2006), a cloud at a distance of 2.5 kpc in the direction towards 
W51, th at was previously stu died in several molecular transi- 
tions by Shirjman et al. (2003), in the 450 //m and 850 //m con- 
tinuum by lOrmel et al.l (120051) . and in the mid-infrared (MIR), 
using the Spit zer Space Telescope, by Kraem er et al. (pro- 
posal 3121, see lVan der Wiel & ShipmanL 12008). These observa- 
tions have revealed several submm condensations and IR point 
sources indicating deeply embedded protostars. 

This configuration is an ideal test bed to distinguish the role 
of internal and external UV fields for the distribution of atomic 
carbon, because the protostars are still too deeply embedded to 
allow their UV radiation to illuminate remote parts of the cloud. 
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We find a clear distinction between inner and outer surfaces. By 
observing the spatial distribution of the [C i] emission and corre- 
lating this with the distribution of molecular material on the one 
hand and with the distribution of young stars within the IRDC 
on the other hand, we can discriminate between the effects of in- 
ternal and external radiation for the formation of atomic carbon. 

In Sect. 2 we describe the observations and the supplemen- 
tary data used in the analysis. Section 3 summarises the results 
and in Sect. 4, we derive abundance ratios for a comparison 
with reference data and give a first interpretation of the results. 
Section 5 summarises the conclusions of the paper. 

2. Observations 

We have used the Sub-Millimeter Array Receiver for 
Two frequencies (SMART, iGraf et al.L 120031) at the Kolner 
Observatorium fur Su b-Millimeter Astronomie (KOSMA, 
IWinnewisser et all 1 19861) 3 m telescope on Gornergrat, Switzer- 
land, to observe the G48 IRDC in the 492 GHz [C i] 3 Pi - 3 P Q 
line. The observations were taken in March 2004 and February 
2005. At that time, the SMART receiver was equipped with 4 
pixels for the 460-490 GHz atmospheric window and 4 pixels 
for the 810 GHz window. The double sideband system temper- 
ature at 492 GHz was about 250 K for three of the pixels, the 
fourth pixel was under-pumped, so that we ignore all data from 
that pixel here. 

We have mapped the cloud in a double beam-switch (DBS) 
raster mode using a chopper throw of 6 arcminutes. Although 
the [Ci] emission is in general extended, we can exclude any 
significant self-chopping effects due to the spatial and velocity 
structure of the IRDC as shown in the Appendix. The HPBW of 
the KOSMA telescope at 492 GHz is 55" and the main beam 
efficiency rj m b =62%. The pointing accuracy of the telescope 
is better than 10", we typically measured a 1 <x pointing scatter 
of about 3" in pointing checks at various sources. The SMART 
pixel separation is 1 15" and we have added in the mapping two 
telescope pointings to fill the pixel spacing so that we obtained a 
slightly undersampled map with a spacing of 38". 

The IF signal was analysed with 4x1 GHz acousto-optical 
spectrometers with a spectral resolution of 1.5 MHz, corre- 
sponding to a velocity resolution of 0.91 km/s at 492 GHz. 
The channel spacing is 1.04 MHz, corresponding to 0.63 km/s. 
The high-frequency channel of SMART was tuned to the upper 
(809 GHz) [Ci] 3 P2- 3 Pi line, but due to the higher receiver noise 
(r sys (DSB) * 650 K) and the lower atmospheric transparency at 
that frequency this provided only an upper limit of about 5 K for 
the upper [C i] line. 

For the analysis, we compare the new [Ci] data data with 
the 13 CO 1-0 maps taken in the Boston University - Five 
College Radio Astronomy Observ atory (BU-FCRAO) Galactic 
Ring Survey (Jack son et al.L 12006). They are given at a velocity 
resolution of 0.21 km/s. While regridding the FCRAO data from 
Galactic to equatorial coordinates, in which the [Ci] map was 
taken, the data were resampled to the same spatial resolution of 
55". We provide all spectra on the original antenna temperature 
scale. 

As both [C i] and 13 CO show emission extending beyond the 
main beams of the telescopes, the measured intensity is given 
by the superposition of the main beam detection and some error 
beam pickup. In the Appendix we have estimated the effective 
beam efficiencies for this configuration using the known main 
and error beam efficiencies and the approximate configuration 
of the IRDC. To convert the antenna temperature to this effec- 
tive beam temperature scale, the KOSMA data have to be di- 
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Fig. 2. Scatter plot of the measured 13 CO and [Ci] intensities 
for different channels and points across the cloud. The red line 
indicates a [Ci]/ 13 CO intensity ratio of 1/3. 

vided by an effective beam efficiency of 0.69 and the FCRAO 
data by 0.70. Since these numbers are almost identical, all ratios 
are unaffected by the beam coupling, so that they have the same 
value when taken on the antenna temperature or the effective 
beam temperature scale. 

3. Results 

Figure [TJ shows all measured [C i] spectra on top of the corre- 
sponding FCRAO 13 CO 1-0 spectra covering the G48 IRDC. We 
find the same spectral features in 13 CO and [C r]. There are four 
dominant velocity components, at 34 km/s, 40 km/s, 52 km/s, 
and 57-59 km/s, but also some emission at intermediate veloci- 
ties between the components at 34 and 40 km/s and the compo- 
nents at 52 and 58 km/s, indicating bridges of gas between the 
components. Within the noise limits we find no significant dif- 
ferences in the spe ctral shapes of the [Ci] and 13 CO 1-0 lines. 
Simon et al. (2006) identified the 34 km/s component with the 
IRDC based on the morphological similarity of the I3 CO compo- 
nent and the mid-IR dark cloud in the MSX image. The 52 km/s 
component has a kinematic distance of more than 4 kpc repre- 
senting a separate cloud which appears accidentally in the same 
line-of-sight but is unrelated to the IRDC. 

Figure|2]shows a scatter plot of the measured [Ci] and 13 CO 
1-0 intensities. The typical [Ci]/ 13 CO intensity ratio is 1/3, the 
vast majority of the points can be described by ratios between 0.2 
and 0.6. We find no systematic deviations of the ratio as function 
of channel velocity, i.e. the same ratio holds in the line centre and 
in the line wings and for most of the velocity components. The 
component at 52 km/s shows some more values at the upper end 
of the range while the component at 58 km/s tends towards the 
lower ratios. Both components are not related to the IRDC as 
they fall at different kinematic distances. At the velocities of the 
IRDC, we find relatively small variations in the ratio between 
the [C i] and the 13 CO intensities for all velocities and all posi- 
tions where the lines are detected. [C i] and 13 CO show the same 
spatial extent, but there are some significant differences in the 
morphology of the emission. In particular, there is remarkably 
strong [C i] emission in the southern part of the cloud, where the 
I3 CO emission is weaker than in the cloud centre. 

In Fig. [3] we show an overlay of the integrated intensities of 
the 34 km/s line in 13 CO and [C i] on top of the recent Spitzer im- 
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13 C0 1-0, CI 3 P,- 3 P , v = 20...70km/s, T A *= -0.2...3.4K 
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Fig. 1. Line profiles measured on 13 CO 1-0 (black, BU-FCRAO Galactic Ring Survey) and in [C i] 3 Pi - 3 P (red, KOSMA) towards 
the IRDC G48.65-0.29. All positions are given relative to the zero position of 19 /! 21 m 44.7", 13°49'31.7" (J2000). For each point 
we show the spectrum covering the velocity range from 20 to 70 km s _I and an antenna temperature range from -0.2 to 3.4 K.The 
positions without [C i] spectra were either not covered in the mapping or covered by a bad pixel providing a noise level much higher 
than the line intensity. The dotted lines at 34 km s _1 mark the LSR velocity of the IRDC. 
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Fig. 3. Integrated intensity contours in I3 CO 1-0 (green) and in [Ci] 3 Pi - 3 fo (yellow) on top of a Spitzer IRAC false-colour 
image, composed of logarithmic intensities in the 8ytzm channel (red), the 5.8//m channel (green), and the 3.6/im channel (blue). The 
contours are labelled with intensitie s [Kkm s" 1 ! inte grated from 30 to 41 km s _1 . The positions of the submm cores PI, P2, and 
EP (from west to east) identified by Orme fetal] (120051) and the southern 24 /zm source S2 from IVan der Wi el & Shipmanl (120081) 
are indicated by asterisks. The circle indicates the beam size of the [C i] observations, the 13 CO data are smoothed to the same 
resolution. 



age of the cloud (IVan der Wiel & Shipmanl 12008). One can eas- 
ily recognise the good correlation between the MIR extinction 
and the 13 CO emission when taking the different angular reso- 
lution into account. The spatial distribution of the [C i] emission 
deviates from this structure, showing two intensity peaks with 
one of them located south of the main IRDC structure at the 
map boundary. 



4. Interpretation 



A full analysis of all the spatial and velocity structure visible 
in Fig. [1] goes beyond the scope of this paper. We provide only 
a first analysis, restricted to the 34 km/s component associated 
with the IRDC and start with some morphological considera- 
tions. 
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Table 1. 

cloud 



Parameters derived at four example positions in the 



Points" 






EP 


PI 


S2 


average* 


Offsets 






0',0' 


1.3',0' 


1.3',-1.3' 




J r b ([Ci])dv 


(Kkms- 


') 


3.1 


4.4 


5.0 


rr 


/ r b ( 13 CO)dv 


(K km s" 


') 


10.4 


11.7 


9.3 


5.1 


r[Q]/r[ 13 co] 






0.29 


0.37 


0.53 


0.31 


N[Ci] c 


(10 lb cm 


2 ) 


5.8 


8.4 


9.4 


3.0 


A?[ 13 CO] 


(10 16 em- 


- 2 ) 


1.2 


1.4 


1.1 


0.6 


N[CO] 


tio" cm- 


- 2 ) 


7.9 


8.8 


7.0 


3.9 


N[H 2 ] 


flO 22 cm 


- 2 ) 


1.0 


1.1 


0.9 


0.5 


N[Ci]/N[CO] 






0.07 


0.09 


0.13 


0.08 



Notes. (n) Within the beam of the [C i] observations, the emission from 
core P2 overlaps with that from PI and EP so that it is not separately 
listed. <6) Spatial average over the [Ci] map. (c) Using the geometric 
mean of emissivities at T„ = 10 K and 25 K. 



lOrmel et~al1 d2005l) identified three submm cores in the 
IRDC as possible sites of early massive star formation. They 
are located at offset positions of 1.26',-0.05' (PI), 0.80',-0.18' 
(P2 ) and Q.05',-0.13' (EP) marked by asterisks in Fig. [3] (see 
also lVan der Wiel & Shin man. 2008). These cores basically line- 
up along the horizonta l filament in Fig. [3] and are also traced 
by the 13 CO emission. IVan der Wiel & Shinmanl (l2008h identi- 
fied 24 yum point sources in all three submm cores (PI contain- 
ing two sources), but also an additional strong 24 fim source 
at 1.35', -1.21' (their S2, visible as blue point source in Fig. [3j, 
south of PI, a position where no submm core was seen. This 
southern part of the cloud, traced by the IR extinction and the 
13 CO emission, is visible in the 850 fim continuum, but shows 
no pronounced condensations. The point sources in the western 
cores (P2 and EP) have not been detected in the IRAC bands 
while the point sources in the eastern core (PI) and the south- 
ern part show a significant flux even at 3.6 /vm with the southern 
source also matching a 2MASS point source. This indicates that 
the latter are either in a more evolved state or at a somewhat 
lower column density relative to the observer. 

Comparing this complementary information with the [Ci] 
emission shows that [C i] neither traces the column density of the 
cloud nor is concentrated at the surfaces of the cloud, but that it is 
peaked at the sites of intense NIR radiation. If we assume that the 
large-scale structure of the IRDC, in particular the surface to vol- 
ume ratio, is similar in all parts of the cloud, the internal, more 
evolved radiation sources remain as the main cause for the en- 
hancement of atomic carbon. The very deeply embedded sources 
present in the western cores do not contribute to a local enhance- 
ment of atomic carbon, but as soon as the sources become visi- 
ble in the NIR, they change the C i abundance. Nevertheless, we 
have to emphasise that we observe relatively small changes of 
less than a factor two in the [C i] intensity. The global C i dis- 
tribution is still dominated by the overall structure of the IRDC 
reflecting the external illumination. 

We can use the available complementary information on the 
structure of the IRDC and on the physics of PDRs to get a rough 
estimate for the temperature of the cloud. In spite of some dif- 
ferences, all PDR models predict that the bulk of atomic car- 
bon is found at temperatures between about 10 K a nd 25 K for 
sources with low or intermediate radiation fields dLe BourlolL 



1 1993b iRollig et"ail l2007lPl Using these values we can deduce 
the C i column density from a single line applying an LTE anal- 
ysis. The LTE emissivity of the [Ci] 3 Pi - 3 Pq transition varies 
between 3.7 x 10~ 17 and 7.6 x 10~ 17 K km/s cm 2 for tempera- 
tures between 10 and 25 K(6.2x 10~ 17 K km/s cm 2 at 15 K) and 
the critical density is about 10 3 cirT 3 . To minimise the relative 
error of the column density due to the temperature uncertainty 
we use the geometric mean of the two extremes, i.e. 5.3 x 10~ 17 
K km/s cm 2 , and obtain a column density estimate that should 
■be accurate within a factor 1 .4. 

TableQ]summarises the results for four sets of effective beam 
temperatures, measured at the locations of the western core (EP), 
the eastern core (PI), the position of the southern 24 fim source 
(S2), and the map average. The observed range of [Ci] bright- 
"ness temperatures is reflected by atomic carbon column densities 
of 3 - 9 x 10 16 cirT 2 . The corresponding beam-integrated optical 
depths are at most 0.25, so that we can neglect self-absorption 
effects. 

Using a similar approach for the 13 CO column densities is 
somewhat less accurate because the molecule exists over a much 
larger temperature range. Here, we can rely on the modelling of 
the G48.65-0.29 IRDC bv lOrmel etail (12001 . They found that 
only a very small part of the submm cores is heated to tempera- 
tures above 20 K while the majority of the gas is subject to tem- 
peratures between 10 and 15 K, i.e. in a range at the lower end 
of the typical temperatures of atomic carbon. This is somewhat 
counter-intuitive, as one would usually expect 13 CO deeper in 
the cloud than atomic carbon at locations better shielded from 
the UV radiation and consequently also significantly colder. 
However, higher densities in the inner cloud also provide a more 
efficient gas-grain interaction that can heat the gas again, even al- 
lowing for non-monotonous temperature distributions. The LTE 
emissivity of the 13 CO 1-0 transition is 7.4xlO~ 16 K km/s cm 2 at 
20 K and 1 .Ox 10~ 15 K km/s cm 2 at 10 K (8.7x 10" 15 K km/s cm 2 
at 15 K), the critical density is well below 10 3 cirT 3 . Using again 
a mean conversion factor, i.e. 8.6 x 10~ 16 K km/s cm 2 should 
be accurat e within 20 %. Assuming a 12 CO/ 13 CO abundance ra- 
tio of 65 (Langer & Penziasl 1 19901) . we can translate this into 
CO column densities and C i/CO abundance ratios. Ap plying the 
standa rd CO/H2 abundance ratio of 8 x 10" 5 from iBlake et al.l 
(fl987l) gives also a rough estimate for the total colu mn densities. 
They a re c onsistent with the n umbers obtained by lOrmel et~afl 
(2005]) and lSimon et al. (2006). All numbers are summarised in 
Table[Q 

We find Ci/CO abundance ratios between 0.07 and 0.13. 
The Ci/CO abundance ratio previously measured in star- 
forming regions covers the range between about 0.03 an d 30 
dMookeriea et alll2006l:lLittle et allll994l:LSun et all l20081. The 
values found for G48.65-0.29 are found at the lower end of this 
range and this is consistent with the picture of atomic carbon 
being produced by star-formation processes in which more qui- 



Some uncertainties exist with respect to the temperature of the 
hottest few percent of atomic carbon, in particular in regions with high 
radiation fields. There, large differences exist between different PDR 
model predictions and this is also the range where the steady-state as- 
sumption used by most models is least justified. A parameter scan using 
the KOSMA-t PDR model showed that higher average temperatures of 
atomic carbon can occur for n/M xx ~ 10 6 ctiT 3 / Mqxo, i.e. for small 
dense clumps in a radiation field exceeding the average interstellar radi- 
ation field xo- The known properties of the IRDC cores are sufficiently 
far apart from that boundary so that the uncertainties do not affect our 
estimate. 

2 Lower values h ave been reported only towards the DR15 Hn region 
bv lOkaetal.l i l2001h . 
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Fig. 4. [C i] 3 Pi - 3 P / ll CO 1-0 line intensity ratio predicted by 
the KOSMA-t PDR model for a spherical, centrally-condensed, 
externally illuminated cloud core of 100 Mq as a function of the 
surface density and the strength of the UV field. Contours are 
drawn at ratios of 0.2, 0.32, 0.5, 1.0, and 3.2 to emphasise the 
region of interest. 



escent clouds, without internal UV fields, show lower Ci/CO 
ratios than more evolved, active star-forming sites. 

The LTE analysis applied above ignores, however, the pos- 
sible spatial separation of the regions dominating the [C i] and 
the 13 CO 1-0 emission. A physically self-consistent model of 
the cloud can be obtained usi ng the known dens ity structure of 
the submm cores obtained by lOrmel et alJ d2005h . The roughly 
spherical cor e structure is well a pproximated by the KOSMA-r 
PDR model (Rollig e t"ai],l2006l) . This model simulates a spher- 
ical cloud with a radial density profile given by a r~' 5 decay 
covering a factor five in radii around a core with constant den- 
sity. Assuming external illumination by an isotropic UV radia- 
tion field and cosmic rays, the model computes the stationary 
chemical and temperature structure by solving the coupled de- 
tailed balance of heating, line and continuum cooling and the 
chemical network us ing the UMIST data base of reaction rates 
dWoodall et all | 2007l) expanded b y separate entries for the 13 C 
chemistry (see Rollig et al., 2007, for details). The model con- 
figuration provides a good match to the quiescent cores in the 
IRDC and may reflect the situation for most of the material in the 
IRDC. Here, atomic carbon is produced in a layer at the surface 
of the dense material and in the more extended diffuse material. 
In the inner core the material is predominantly molecular. 

This approach is, however, not appropriate for the cores with 
internal radiation sources, PI and S2, where we expect some ad- 
ditional atomic carbon production in the centres. Consequently, 
the KOSMA-r PDR model should be able to reproduce the [C i] 
3 fi - 3 P / 13 CO 1-0 intensity ratio of about 0.3 found at the ex- 
tinction peak and in the cloud average, but will only provide a 
lower limit to the [C i]/ 13 CO ratio for the two warm cores. When 
comparing the line ratios at the different positions (see Tab. Q]), 
we find that the internal radiation from PI produces only a slight 
enhancement of the atomic carbon abundance over the value pro- 
duced by the external UV illumination while for the southern 
source internal and external illumination provide similar contri- 
butions. 



Figure |4] shows the predicted intensity ratios as a function of 
cloud surface density and radiation field provided by the PDR 
extraction too0 for the known core mass of about 100 Mq 
and solar metallicityQ. White contours indicate the area of ra- 
tios matching the observed range. The ratio of about 0.3, i.e. 
lg{ J T([Cl])dv/ J T( 13 CO)dv] * -0.5, observed within most of 
the cloud and marked by the central white contour is only ob- 
served at high densities. For radiation fields in the order of th e 
standard interstellar radiation field xo (Draine & BertoldiifT9%h . 
a model with a clump surface density of about 3 x 10 4 cirT 3 re- 
produces the observed ratio. The model predicts average temper- 
atures of atomic carbon and 13 CO of 12 K and 9 K, respectively, 
consistent with the estimates used above for the LTE analysis, 
but somewhat at the low edge as the model does not contain in- 
ternal UV radiation sources. 

When comparing the model with the r esults from the fit o f 
the dust continuum and HCO + emission bv lOrmel et all (120051) 
we fin d a major discrepancy in the gas density. Orme Letdl 
(2005) derived surface densities of only about 10 3 cm -3 . When 
taking internal dumpiness of the cores into account, local den- 
sities grow by a factor ten (core PI) or five (P2 and EP), but 
are still significantly lower than the model density fitting the 
[C i]/ 13 CO ratio. For the lower densities, the KOSMA-t model 
always predicts significantly higher line ratios. Models that as- 
sume a further break-up of the cores into smaller clumps have 
a higher surface to volume ratio, thus predicting even higher 
[Ci]/ I3 CO ratios. Lower ratios could be achieved for radiation 
fields significantly below the standard field Xo but that would be 
hard to explain. Excitatio n from W5 1 can be excluded due to 
the large relative distance dReTdl[T99l so that we expect fields 
in the order of xo- A match could also be achieved when as- 
suming unrealistically high metallicities Z > 3. With realistic 
parameters and the known masses and densities, the model pre- 
dicts ratios above 0.5, only observed at the peaks with internal 
radiation sources where the model is expected to underestimate 
the ratio by neglecting this contribution. Invoking other model 
geometries than spherical clumps can only worsen the situation 
as they will always increase the surface to volume ratio of the 
cloud. Fractionation effects will hardly resolve the discrepancy 
as the KOSMA-t PDR model includes the full 13 C chemistry. 
But we have to admit that reaction rates for isotope-selective re- 
actions are quite uncertain at low temperatures (see Rollig et al., 
in preparation). Thus, there is a clear contradiction between the 
measured [C i]/ I3 CO line ratios and the model predicting higher 
values. 

The same contradiction was noticed bv lBensch et all (|2003) 
for the MCLD 123.5+24.9 translucent cloud. They suggested a 
pre-shielding by dust in a thin ionised gas as a possible solution. 
Time-dependent chemistry has been shown to ten d towards the 
opposite direction, incr easing the C i abundance dStorzer et all 
1997; lOka et all 120041) . but a fully time-dependent model, in- 
cluding turbulent mixing, still awaits development. 

5. Conclusions 

Our observations of [Ci] in an IRDC showed an overall very 
good correlation between atomic carbon and molecular gas, seen 

3 http://hera.phl .uni-koeln.de/~pdr/ 

4 The web tool provides integrated line intensities in terms of photon 
energies, i.e. units of erg s -1 cirT 2 sr~'. For intensity units of K km s -1 , 
the conversion factor of 2fcv 3 /c 3 = 1.025 x 10~ 15 (v/GHz) 3 has to be ap- 
plied. This gives a conversion factor for the j T([Cl])dv/ j T( l3 CO)dv 
intensity ratio of 0.01 12. 
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in 13 CO 1-0, both in the spatial and the velocity structure. This 
is naturally explained by a fracta l cloud structure illumi nated 
by the external Galactic UV field dOssenkopf et al.L 12007). The 
quantitative [Ci] 3 P\ — 3 Po/ 13 CO 1-0 ratios, however, cannot 
be explained by available steady-state models. A spherical PDR 
clump model predicts higher line ratios than observed. In a frac- 
tal structure with more internal surfaces, this ratio should even 
grow further. More sophisticated models are needed. Future ob- 
servations of the emission of atomic oxygen at 4.7 THz may help 
to resolve the puzzle as [O i] is another, independent tracer of il- 
luminated, dense cloud surfaces. 

With the significant enhancement of the [C i] intensity ob- 
served towards the southern source, visible in the Spitzer IRAC 
bands, we can provide a first distinction between the roles of 
inner and outer PDRs if we assume that the outer structure is 
similar in all parts of the IRDC. A fraction of the [Ci] emission 
comparable to that from the outer cloud surfaces, exposed to the 
interstellar radiation field, can be attributed to the internal radi- 
ation source. This means that any young stellar object that has 
created a sufficient cavity in the parental cloud, turning visible 
in the NIR, produces a considerable fraction of the total [Ci] 
emission. For most of the star-forming regions observed so far 
in [Ci], the emission will thus stem from inner surfaces. 

Acknowledgements. We thank an anonymous referee for many helpful com- 
ments. This work has been supported by the Deutsche Forschungsgemeinschaft 
through grant 494B. It has made use of NASA's Astrophysics Data System 
Abstract Service. 



Appendix A: Effects of extended emission 

Both 13 CO and [C i] are known to show very extended emission 
so that it needs to be investigated to which degree the measured 
intensity at each point may be affected by this emission. Two ef- 
fects have to be considered: self-chopping and intensity pick-up 
by the extended error beam. Taken the similarity of the emis- 
sion of 13 CO and [Ci] we can estimate these effects from the 
larger 13 CO 1-0 FCRAO map showing also the environment of 
the IRDC. It is displayed in Fig. lA.ll using the same integration 
range as in Fig. [3] 

The [Ci] observations were taken in double beam-switch 
mode using a fixed chop throw of 6' in azimuthal direction. The 
sky rotation relative to the azimuthal direction then leads to vari- 
able OFF positions depending on date and time of the observa- 
tion. In Fig. I A. 1 1 we have indicated the area of OFF maps cov- 
ered in the course of the observation runs. We see that the centre 
of the map is never affected by self-chopping effects, but that 
the north-west corner and the east edge of the map could be af- 
fected by self chopping. As we have not detected any significant 
[C i] emission at these points, this might only mean that [C i] is 
slightly more extended than shown in Fig. [3] 

We have to admit, however, that the spectra show some self- 
chopping effects at frequencies out of the considered velocity 
interval. This can be seen as negative features around 44 km/s 
in some spectra in the south-west part of the map and around 
55 km/s in the south of the map. These contributions stem from 
gas components south-east of the IRDC. As all these features 
only show up outside of the velocity integration interval assigned 
to the IRDC, self chopping effects can be neglected for the inte- 
grated [Ci] intensities. 

The extended emission will also lead to additional intensity 
pick-up in the error beams of the telescopes. The error beam of 
the KOSMA telescope, produced by a slight misalignment of th e 
individual panels has a size of about 330" (Kramer et al., 2004). 




290.6 290.5 290.4 290.3 290.2 
RA ['] 



Fig.A.l. Map of the large-scale emission as measured by the 
13 CO 1-0 line. The map covered by [Ci] observations is shown 
as gray rectangle in the centre, the brown rectangles around that 
map indicate the locations of the OFF maps taken in the dou- 
ble beam-switch scheme at the different times and dates of the 
observations. The circle indicates the size of the KOSMA error 
beam at the [C i] frequency. The FCRAO error beam at the 13 CO 
1-0 frequency is smaller by a factor three and is not shown here. 



Comparing the moon efficiency of 80 % to the main beam ef- 
ficiency of 62 % provides a maximum error beam efficiency of 
18 %. We have estimated the error beam pickup by convolving 
the FCRAO map with a 330" Gaussian, weighted by the error 
beam efficiency, and comparing it to the map values weighted by 
the main beam efficiency. We find that at the emission peaks, the 
error beam contribution falls between 8 and 12%, over most of 
the map it amounts to about 15 % but it reaches a value of 30 % 
at the south-west corner of the map. Consequently, the true [C i] 
emission might be somewhat less extended in the south-west 
direction than shown in Fig. [3] Due to the error-beam pickup 
we overestimate the local emission when dividing the measured 
intensity by the main-beam efficiency only. Using an effective 
beam efficiency of 69 % instead of 62 % gives a correction for 
this effect that should be accurate to about 10 % over most of the 
observed map. 

The error beam of the FCRAO telescope at the 13 CO 1- 
frequency has a width of about 100". The main beam effi- 
ciency is about 48% and the measured moon efficiency 72 % 
dLadd & Covevi 120001) . Repeating the deconvolution exercise 
for this smaller error beam gives a relatively constant ratio of 
about 0.45 for most of the map which is only slightly smaller 
than the ratio between error beam and main beam efficiency. 
Consequently, we will use a value of 70 %, only slightly below 
the moon efficiency as the effective beam efficiency to calibrate 
the I3 CO. 

This approach may slightly underestimate the intensity as 
the moon efficiency is only an upper limit for the sum of main- 
beam and error-beam efficiency and some small components at 
the scale of a few arcminutes might contribute, but the residual 
error from this contribution will certainly fall below 10 %. 
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